Possible centers of broadband near-IR luminescence 
in bismuth-doped solids: 
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Subvalent bismuth centers (interstitial Bi + ion, Big + cluster ion, and Bi2 cluster) are examined 
as possible centers of broadband near-IR luminescence in bismuth- doped solids on the grounds of 
quantum-chemical modeling and experimental data. 
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I. INTRODUCTION 

Near-IR broadband (1100 - 1400 nm) luminescence 
in bismuth-doped glasses discovered in Ref. [ij is be- 
ing studied intensively. By now the luminescence has 
been observed in many bismuth-doped glasses, such as 
alumosilicate (e.g. [3-Q), alumogermanate (e.g. 0, 0~ 
9|), alumoborate (e.g. @, [13]), alumophosphosilicate, 
alumophosphate, alumophosphoborate (e.g. pj. [Til fl^|). 
chalcoeenid e fi~3l . Q , in several bismuth-doped crystals 
(RbPb 2 Ck [lf7 2MgO • 2A1 2 3 • 5Si0 2 cordierite Jlp , 
BaB 2 4 BaF 2 [ljLBa 2 P 2 7 [ljj, Bi^AlCLOg [If), 
and in FAU-type (e.g. |21|) and Y-type [22| zeolites. The 
bismuth-related IR luminescence is used successfully in 
laser amplification and generation (e.g. dHHU). 

However there is no commonly accepted model of the 
IR luminescence center. Several models are suggested, 
such as electronic transitions in Bi interstitial atoms 
P, in Bi+ [TMESES Bi 2+ [H and Bi 5+ 
interstitial ions, in BiO interstitial molecules [4J1, Bi 2 , 
Bi^ and Bi 2- interstitial dimers 0, 0, H Elllj, in 
other neutral bismuth clusters 0], in Big + cluster ions 
H HE Hi, Bi0 4 complexes with tetrahedral co- 
ordination of the central bismuth ion [29| . in oxygen- 
coordinated complexes formed by pairs of threefold coor- 
dinated bismuth atoms [30] , in complexes formed by Bi + 
or Bi 2+ substitutional or interstitial ions with certain de- 
fect centers, such as anion vacancies [3lj |. 

Recently subvalent bismuth centers such as Bi + single- 
charged bismuth ion and Big + cluster ion have attracted 
considerable attention as possible sources of the near-IR 
luminescence [20, 26, 28]. In this paper we examine three 
subvalent bismuth centers basing on experimental data 
available and our quantum-chemical modeling. Namely, 
the above-mentioned (interstitial) Bi + ion and Big + clus- 
ter ion, and (as a novel suggestion) Bi" cluster are stud- 
ied. 
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II. CALCULATIONS 

Configuration and electronic states of Big + and Bi" 
clusters were calculated using GAMESS (US) quantum- 
chemical code [32]. Firstly, ground state DFT calcu- 
lations of the clusters were performed using B3LYP1 
functional (Becke+Slater+Hartree-Fock exchange and 
Lee-Yang-Parr+Vosko-Wilk-Nusair (5) correlation — see 
GAMESS (US) manuals [H for details) and effective 
core potentials (ECP) and bases of three types devel- 
oped in Refs. 13314351 with one Huzinaga's polarization 
d-type function [36j added in each basis. 

The calculated molecular orbitals were then used to 
obtain the final molecular wave functions by means 
of perturbation theory with multiconfigurational self- 
consistent-field reference functions. These final wave 
functions and corresponding eigenvalues were used to cal- 
culate eigenfunctions and energies of the ground and ex- 
ited states of the clusters by method of configuration in- 
teraction in active orbitals space with spin-orbit interac- 
tion taken into account. 

Multiconfiguration complete active space self- 
consistent field (MC-CASSCF) and multircference 
configuration interaction (MRCI) were used to calculate 
electronic states of Big + and Bi" clusters. Separate 
CASSCF calculations were performed for the ground 
state and the lowest excited state of each cluster. Then 
MRCI calculations were made with all CASSCF configu- 
rations with coefficients large enough taken into account 
(see GAMESS (US) manuals for details). Additionally, 
single and double excitations from these configurations 
were allowed. 

Evaluative calculation of optical spectra in external 
electric field and of ground-state optical absorption spec- 
tra (to verify independently the results of configura- 
tion interaction calculations) were performed by time- 
dependent density functional theory (TDDFT) method 
using Octopus program [371 and Hartwigsen-Goedecker- 
Hutter pseudopotentials [38| with spin polarization and 
spin-orbit interaction taken into account. PBE density 
functional [391 ] was used in the ground-state calculation. 
In the Octopus code, to obtain the linear optical absorp- 
tion spectrum of the system, all frequencies of the sys- 
tem are excited by giving certain (small enough) momen- 
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turn to the electrons. Then the time-dependent Kohn- 
Sham equations evolves in real space for a certain real 
time and the dipole-strength function (or the photo- 
absorption cross section) is obtained by a Fourier trans- 
form of the time-dependent dipole moment. Adiabatic 
LDA approximation is used in these calculations to de- 
scribe exchange-correlation effects. The Octopus code 
uses real-space uniform grid inside the sum of spheres 
around each atom of the system (a single one in our case). 
The sphere radius and the grid spacing were taken to be 
was 8.0 and 0.25 A, respectively, in our calculations. The 
real-time propagation was performed with 2 • 10 4 time 
steps with the total simulation time of about 20 fs. The 
Fourier transform was performed using third-order poly- 
nomial damping (see (37| for details of the code). 



III. ELECTRONIC STATES AND OPTICAL 
SPECTRA OF FREE BISMUTH SPECIES 

A. Bi + bismuth ion 

The optical absorption and emission of free singly 
charged positive bismuth ion, Bi + , are well known (see, 
for example, |4ll|). The ground state and two low- lying 
excited states with energies about 13325 and 17032 cm -1 
are known to be formed owing to spin-orbit splitting from 
the 3 P term. The electronic configurations of these three 
states are 6p^ 2 (J = 0), 6pJ/ 2 +6pg/ 2 (J = 1), and 

6pjy 2 +6pg/ 2 (J = 2), respectively. The subsequent ex- 
cited states with much greater energy (about 33940 and 
44175 cm -1 ) correspond to X D and 1 S terms with elec- 
tronic configurations 6pgy 2 (J = 2) and 6p^ 2 (J = 0), 
respectively. Obviously, all the electric dipole (El) tran- 
sitions between the states just listed are forbidden, only 
magnetic dipole (Ml) and electric quadrupole (E2) tran- 
sitions being allowed. In particular, the ground state 
of Bi + ion is characterized by weak absorption at wave- 
lengths about 750, 585 and 295 nm corresponding to Ml 
and E2 transitions from the ground state to three lowest 
excited states. The most long- wavelength luminescence 
in Bi + ion corresponds to the Ml transition from the first 
excited state to the ground one (Fig. [TJ (a)). Obviously, 
there is no Stokes shift in free ion. 

In Ref. [42| absorption spectra of Bi + ion in chlorides 
melts have been studied and interpreted basing on crystal 
field theory. According to this study, two lowest excited 
states of the ion are split by the crystal field in two and 
three levels, respectively. As a result, the transition en- 
ergy decreases significantly, and bands around 900, 785, 
660-690, and 585 nm are observed in the absorption spec- 
trum (Fig. [TJ (b)). 

Unfortunately, calculation of Bi + ion spectra (both free 
ion and one in electrostatic field) with the spin-orbit in- 
teraction taken into account presents severe computa- 
tional problems using GAMESS (US) program. Our con- 
servative estimates and evaluative calculations performed 
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Figure 1. Experimental states and transitions in (a) free Bi + 
ion [!H and (b) Bi+ ion in AlCl 3 -NaCl eutectic gg]. 



by the TDDFT method using Octopus code causes us 
to anticipate that the wavelength of the lowest-energy 
transition does not exceed 1000 nm essentially. Such a 
conclusion is consistent with the results of Ref. [421 ■ 



B. 



Bi| + cluster ion 



Our calculations prove the electronic ground state of 
Big" 1 " cluster ion to be a singlet. In this state, the ion 
symmetry is found to be D 3?l 
distances Bi ax 



Bi eq and Bi eg 



(Fig. |21 ((a) I). Interatomic 



Bi e „ between axial, Bi Q 



and equatorial, Bi eg , bismuth atoms are 3.089 A and 
3.325 A, respectively. This agrees well with both ex- 
perimental and calculated data for free Big 4 " cluster ion 
[48|] . In Fig. [3] (a) are shown the calculated scheme of the 
Bi 5 + cluster ion levels and transitions between them. In 



3 





(a) 
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Figure 2. Calculated configurations of free Bii? 4 " cluster ion: (a) ground-state singlet (D 3h symmetry), (b) first excited-state 
triplet (C 2v symmetry). 



the ground state configuration, the electronic states are 
classified according to irreducible representations of the 
D 3 ^ group. From the ground state with the A 1 symmetry 
El transitions are allowed to excited states with E' and 
A 2 symmetry, giving rise to absorption at 750, 585 and 
450 nm wavelengths (shown by solid arrows in the Fig. [3]). 
When perturbed by an external electrostatic (crystalline) 
field, El transitions from the ground state to the excited 
states with E symmetry become slightly allowed corre- 
sponding to absorption at wavelengths of about 835 and 
730 nm (displayed as dashed arrows in Fig. [3]). It should 
be noted that in the absence of spin-orbit interaction all 
the excited states in Fig. [3] correspond to the first ex- 
cited triplet. These results are in good agreement both 
with the experimental data available on Big + cluster ions 
optical absorption [48|, [S(J, and with the results of our 
TDDFT calculations using the Octopus code shown in 
Fig. [4] As well those agree reasonably with experimental 
data for Bi 5 (AlCl 4 ) 3 crystal [|o|. 

Fig. [2] ( |(b)[ ) shows the calculated configuration of Big + 
cluster ion in the lowest excited state. To our knowledge, 



the calculation of configuration of this ion in an excited 
state is performed for the first time. Being excited to 
this state, Bii? + cluster ion undergoes deformation cor- 
responding to E" irreducible representation of the D 3fl 
group. As a result, the cluster ion symmetry is reduced 
to C 2 „. One of the Hi ax — Bi eg interatomic distances is in- 
creased to 3.363 A, the other two distances are decreased 
to 3.043 A, and corresponding Bi e(3 — Bi eg distances be- 
come 3.332 A and 3.362 A, respectively. Calculated levels 
of Big + cluster ion in the configuration of its lowest ex- 
cited state, transitions to the ground state and relative 
lifetimes are shown in Fig. [3] (b). In our calculations 
two El transitions to the ground state are found in this 
configuration. Those lead to two bands of IR lumines- 
cence in 1200 - 1290 nm range with lifetimes differing 
by an order of magnitude, which may be called "slow" 
and "fast" components (unfortunately, only relative life- 
times may be estimated in our calculations) . It should be 
remarked that such two components of bismuth-related 
IR luminescence have been observed repeatedly (see e.g. 
Refs. 
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Figure 3. Calculated states, excitation transition rates and 
luminescence lifetimes in free Big + cluster ion: (a) ground 
state, (b) 1-st excited state (transition rates and lifetimes in 
relative units). 



Furthermore, there are luminescence bands near 1000 
and 900 nm (as well "slow" and "fast", respectively). 
These results are also confirmed in our TDDFT calcula- 
tions. 



C. Bi2 cluster 

According to our calculation the electronic ground 
state o f Bi° cluster turns out to be a singlet. As seen from 
Fig. [5] ( (a) ), in this state the Bi° cluster is a regular tetra- 
hedron (T d symmetry). Bi— Bi interatomic distances are 
found to be 3.104 A in agreement with the results of 
Ref. (46|. Fig. |5] (a) presents the calculated scheme of 
levels, transitions between those and the relative inten- 
sity of the transitions in the Bi° cluster. The electronic 
states of the cluster in the ground state configuration cor- 
respond to T d group irreducible representations. Thus, 
the ground state has A 1 symmetry From the ground 
state El transitions only to the excited state with F 2 
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Figure 4. Absorption cross section of (a) Big + and (b) Bi2 
clusters calculated by TDDFT method. 



symmetry are allowed. Those correspond to absorption 
bands near 670 and 420 nm (solid arrows in Fig. \6\) . In 
an external electrostatic field El transitions from ground 
state to the states with A 2 , F l7 and E symmetry become 
slightly allowed as well, and absorption at wavelengths 
of about 735, 590 and 525 nm arises (dashed arrows in 
Fig. IH]). The first four excited states are formed due to 
spin-orbit splitting of the 3 F 2 triplet state, and the next 
ones arise from 1 F 2 and 1 F 1 singlet states. 

The energies of states and transitions in Bi" cluster 
obtained in our calculation agree well with the results of 
the calculations in Ref. [46 : ] . However the results differ 
in the order of levels: the states with E and A 2 symme- 
tries turn out to be interchanged. On the other hand, the 
transition energies obtained in quantum-chemical calcu- 
lations, both our and that of Ref. [46| . agree closely with 
our TDDFT calculation (Fig.gJ. ' 

Excitation to the lowest excited state with E symme- 
try is accompanied by rearrangement of the cluster con- 
figuration resulting in the cluster symmetry reduced to 
C 2l) . The length of one of the tetrahedron edges (Bi— Bi 
distance) is increased to 3.378 A and the length of the op- 
posite edge is reduced to 3.005 A, while the other Bi— Bi 
distances are reduced to 3.082 A (Fig. [5] ([(b)])). 

In Fig. [5] (b) are shown the calculated levels together 
with transitions to the ground state and their relative 
lifetimes for the lowest excited state configuration of Bi" 
cluster. According to the calculations, in this configura- 
tion the most long-wavelength luminescence in the 950 
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- 1040 nm range contains two components ("slow" and 
"fast" in the same degree as above) corresponding to El 
transitions from the first and second excited states to 
the ground one. Furthermore, El transitions from two 
higher excited states to the ground state can lead to the 
"fast" luminescence near 830 and 660 nm. Our TDDFT 
calculations with the Octopus code yield similar results. 



IV. BISMUTH SPECIES IN SILICATE GLASS 
NETWORK 



To study the possibilities of occurrence of Bi ion, Big 
cluster ion, and Bi^ cluster in silicate glass network, we 
performed quantum-chemical modeling of these centers 
in ring interstitial sites of silica glass network using clus- 
ter approach. In the clusters there are either one or two 
six-member rings formed by four Si04 tetrahedra. Dan- 
gling bonds of the outer oxygen atoms in the clusters were 
saturated with hydrogen atoms. The bismuth species 
was placed initially at the center of the ring(s) and then 
the complete geometry optimization was performed. De- 
tails of computational techniques were described above. 
Standard 3-21G basis was used for saturating hydrogen 
atoms. 



A. Bi bismuth ion 

In our calculations equilibrium position of Bi + ion in 
silicate glass network is found in interstitial sites formed 
by six-member rings, Bi + ion being located near the rings 
axis between the rings planes. However, the calculations 
prove this position of Bi + ion to be not stable enough: 
even relatively small displacement of the ion (not more 
than 10 % of the minimal distance to neighboring atoms) 
leads to further drift apart from the equilibrium position. 
Subsequent process is determined by the environment. 
Thus, in the regular ring interstitial Bi + ion moves to- 
wards one of the nearest oxygen atoms in the surround- 
ing rings and threefold coordinated bismuth atom and 
threefold coordinated silicon atom are formed. If there 
are another Bi + ion(s) or Bi° atom(s) nearby, interstitial 
dimers or larger bismuth clusters are likely to arise in- 
stead. As an example, calculated position of BiJ dimcr 
in six-member rings interstitial of alumosilicate network 
is shown in Fig. [7J Some evidence is found in our calcu- 
lations that Bi + ion may be stabilized in interstitial site 
of the silicate glass network by neighboring defects, such 
as oxygen vacancies. 

Thus, although single interstitial Bi + ions may occur 
in the silicate glass network, their stability is poor and 
significant concentration of such ions seems hardly prob- 
able. On the other hand, small bismuth clusters or Bi + 
— defect complexes may be formed in the network intcr- 
stitials, such clusters and complexes being quite stable. 

It is notable that the luminescence Stokes shift which 
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Figure 6. Calculated states, excitation transition rates and 
luminescence lifetimes in free Bi2 cluster: (a) ground state, 
(b) 1-st excited state (transition rates and lifetimes in relative 
units) . 



vanishes in free Bi + ion should be expected to be low in 
the interstitial ion as well, but this is not the case for the 
above-mentioned clusters or complexes. 



B. Bis+ cluster ion 

Geometric sizes of Big + cluster ion are obviously too 
large to allow the ion to be incorporated in ring intersti- 
tial sites of silicate glasses. On the other hand, owing to 
significant electric charge the ion is very likely to react 
with the neighboring electronegative atoms. 

Our calculations performed using the above-described 
cluster models confirm these considerations. It is found 



5 cluster ion reacts with the bridging oxygen 



that Bi 

atoms of the network to form interstitial B14 cluster and 
single Bi 3+ ion nearby. The latter is then built into 
the glass network forming threefold coordinated bismuth 
atom. The resulting structure is shown in Fig. [8l 

It should be realized that even Bii? + ion is not able to 
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Figure 7. Calculated configuration of Bi^~ dimer in six- 
member ring interstitial of alumosilicate glass network. 



react with neighboring network atoms for whatever rea- 
son, the 3+ charge state of the Bis cluster turns out 
to be unstable with respect to electron capture. We 
have calculated electron affinity of Big + cluster ion by 
the above-described quantum-chemical methods in DFT 
level of theory and obtained closely agreed values for ver- 
tical and adiabatic electron affinities: 14.2 eV < EA V < 
EA a < 14.6 eV. 

These values significantly exceed not only the elec- 
tron affinity of the silicate network (< 1.2 eV), but its 
band gap width (< 9 eV) as well. Hence Big + cluster 
ion should capture electron in silicate host, the capture 
being possible not only from the conduction band, but 
from the valence band as well. 

Thus, interstitial Big + cluster ions occurrence should 
be recognized to be impossible in silicate hosts. 

However, it seems reasonable to guess that Bi| + clus- 
ter ion can occur in zeolite hosts due to extremely large 
size of their cage-like interstitials and stabilizing action of 
(AIO4) 1 tetrahedra similar to such action of (AICI4) 1- 
ions in Bi5(AlCi4) 3 crystal [20l |. 



C. Bi2 cluster 

The calculations have proved B14 clusters to be able 
to occur in ring interstitials of silicate networks, despite 
the large size of the cluster. One such example is above- 
described formation of B14 cluster from Big + cluster ion. 
However threefold coordinated bismuth atom does not 
need to be present in the surrounding rings: according 
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Figure 8. Calculated configuration of Bi^ cluster in six- 
member ring interstitial of silica glass network. 

to our modeling, a configuration quite similar to that 
shown in Fig. [8] arises as well in purely-silicate and alu- 
mosilicate rings. Such positions of B14 clusters in ring 
interstitials are found to be extremely stable. Even rel- 
atively large displacement and (or) rotation of the Bi° 
cluster, accompanied by its atoms shifted up to 50 % of 
interatomic distance do not make the above-mentioned 
configuration be destroyed. 

Thus, B14 clusters can occur in the ring interstitial of 
the pure-silicate network, and probably in certain binary 
silicate networks as well, such an interstitial position be- 
ing extremely stable. Together with the fact the Bi° clus- 
ter binding energy is known to be abnormally high [47| . 
this suggests the possibility of high concentration of such 
clusters in silicate hosts. 



V. CONCLUSIONS 

Both the experimental data available and the results 
of our calculations allow draw some conclusions concern- 
ing the origin of the above-mentioned IR luminescence 
centers in the bismuth-containing systems. 

1. For Bi + singly-charged positive ion: 

(a) the lowest-energy transition wavelength does 
not exceed 1000 nm; 

(b) all the relatively long-wavelength transitions 
occur between electronic states originating 
(due to spin-orbit interaction) from the same 



3 P term and hence intensities of both absorp- 
tion and luminescence transitions are of the 
same order; 

(c) Stokes shift is low or even vanishes; 

(d) although there are equilibrium interstitial po- 
sitions of Bi + ion in certain solids, single Bi + 
seems to be not stable enough, at least in ox- 
ide glasses networks. 

Thus, single Bi + ion can be responsible for certain 
optical absorption bands in bismuth-containing 
systems, but can not be considered as a center of 
the IR luminescence. 

2. For Big + cluster ion: 

(a) IR luminescence spectrum accords well 
enough with the experimentally observed lu- 
minescence in bismuth-containing systems; 

(b) there are a number of bands in the absorption 
spectrum corresponding to the absorption ob- 
served in experiment; 

(c) in the absorption spectrum of free Big + clus- 
ter ion there is no band at 500 nm, the most 
intense in experimental spectra of bismuth- 
containing systems; 

(d) large size, high electrical charge and high elec- 
tron affinity of the Big + cluster ion make it 
practically impossible for the ion to occur in 
a solid host other than, supposedly, zeolite. 

(e) in silicate glass the cluster ion can react with 
bridging oxygen atoms of the glass network 
forming threefold coordinated bismuth atom 
bound in the network and interstitial B14 clus- 
ter. 

Thus, Big + cluster ion can be responsible for a num- 
ber of absorption and IR luminescence bands in 
certain bismuth-containing systems, namely, ionic 
solutions and, possibly, zeolites, but not for those 
in silicate glasses. 

3. For Bi" cluster: 

(a) in free Bi° cluster, IR luminescence occur only 
near 1000 nm wavelength, but under the in- 
fluence of crystal field the wavelength can in- 
crease (up to 1100 - 1150 nm according to our 
estimations); 

(b) in the absorption spectrum of the cluster there 
are a number of bands corresponding to the 
absorption observed in experiment including 
the band at 500 nm. However, the band near 
670 nm is the most intense, which is hardly 
consistent with experimental data; 

(c) there are equilibrium positions of B14 cluster 
in sufficiently large interstitial sites in solid 
hosts (e.g. in the ring interstitials of oxide 
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glasses networks), these positions being quite 
stable. 

Thus, Bi° cluster can be responsible for a num- 
ber of absorption and IR luminescence bands in 
the bismuth-containing systems. 

4. If in a bismuth-containing system both Bi + ions 
and Big + cluster ions or (and) Bi° clusters occur 
simultaneously, one would expect mutual excitation 
transfer between the centers of different types. 

Thus, none of the Bi + , Big + , and Bi° centers alone can 
not give rise to the experimental IR luminescence and 
absorption spectra typical for the bismuth-doped solids. 
On the other hand, these spectra of ionic solutions are 
likely to be caused by Bi + and Big + ions together. 

It is our opinion that the experimental data available 
are consistent with the following speculations involving 
the centers under consideration. In solid hosts, most 
likely both Bi + ions and Bi° clusters occur in interstitial 
sites (the latters in sufficiently large interstitials, in the 
ring ones particularly). Moreover, it is conceivable that 



interstitial Bi + ions could form pairs with each other or 
with interstitial Bi° atoms. In ionic solutions, Bi + ions 
and Big + cluster ions most likely occur. Excitation trans- 
fer from Bi + ions to other bismuth centers is expected to 
proceed in all the cases. 

On the other hand, it should be emphasized that no 
experimental data nor calculation results suggest that in 
the bismuth-doped systems there are no IR luminescence 
centers other than those considered in this paper. In par- 
ticular, we persist in the belief that negatively charged 
bismuth dimers, Bi^~ and Bi^ - [l6L [27l|. and dimers (or 
even more complex clusters) of Bi + or Bi 2+ ions and 
some intrinsic defect centers [3l[ are quite possible to 
occur in bismuth-doped solids. Moreover, it seems rea- 
sonable to expect that bismuth centers of several types 
may be formed concurrently in these solids. The assump- 
tions concerning existence of several types of IR lumines- 
cence centers in bismuth-doped solids is supported by 
significant variation of both absorption spectra and IR 
luminescence spectral and temporal dependencies, with 
preparation conditions etc. observed in various systems. 
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